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Abstract 

We present a simple but general treatment of neutrino oscillations in the framework of 
quantum mechanics using plane waves and intuitive wave packet principles when necessary. 
We attempt to clarify some confusing statements that have recently appeared in the literature. 
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The quantum mechanics of neutrino oscillations have been studied in several papers 
0, 0, |, ^, H 1^, H, 1^, [ly, [Tl[] and reviewed in [0, . Alternative derivations of neutrino 
oscillations in the framework of quantum field theory have also been presented in |15|, 
161 , p!7| , |T8|, |1^, |20|, |21|, |2^, Q . In spite of all these studies, the understanding of the theory 
of neutrino oscillations seems still unsettled, as one can see from controversial claims in the 

M m m 



recent literature El, |25 



In this article we wish to clarify some of these issues using a simple quantum mechanical 
treatment of neutrinos as plane waves. A rigorous treatment of neutrino oscillations in the 
framework of quantum mechanics requires a wave packet formalism 0. Here we employ 
some of its principles only when necessary, without technical details (see [Q, ^ |T^). 

Neutrino oscillations are a consequence of neutrino mixing and the fact that neutrino 
masses are very small. Neutrino mixing implies that the left-handed components u^l of the 
flavor neutrino fields (a = e, /i, r)[] are superpositions of the left-handed components u^l of 



neutrino fields with definite mass {k 



1,2, 



3): 



k=l 



ak J^kL 



(a 



e,iJ,,T) 



^ We consider for simplicity only the three active neutrino flavors whose existence is firmly established (see, 
for example, Ref. [|3l|]). However, the formalism is valid for any number of neutrinos, including additional 
sterile neutrinos whose existence is under investigation (see, for example, Refs. js^ , |3^). 

The flavor of an active neutrino {ve, t'^, i't) is determined by the associated lepton (e, fi, t) in charged- 
current weak interactions. Sterile neutrinos do not feel weak interactions, as well as electromagnetic and 
strong interactions (as active neutrinos); they are sensitive only to gravitational interactions. 
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where U is a unitary matrix (f/f/'^ = f/'^f/ = 1), i.e. such that 

'^UakUp,. = 6ai3 , (2) 
k 

UakUaj = 5kj ■ (3) 

a 

A neutrino with definite flavor is described by the state 

Wa) = Y.u:kWk), (4) 

k 

such that 

(0| Ua = ^ UakUpj (0| Uk Wj) OC ^ UakU^k = , (5) 

k,j k 

where we used (0| {I'j) oc Skj. The state describing a flavor antineutrino has complex- 
conjugated elements of the mixing matrix with respect to the flavor neutrino state in Eq. (^): 

\i^a) = ^Uak \ i^k) ■ (6) 

k 

This implies that the oscillations of neutrinos and antineutrinos are different only if the 
mixing matrix is complex {i.e. only if there is CP violation). 

Let us consider z/„ produced at time tp and space coordinate^ xp by a weak interaction 
process with a speciflc flavor a {a = e, fi, r). This neutrino is described by a state \ip^) that 
is a superposition of neutrino states with deflnite kinematical properties: 

\£) = Y.^:kWk,Pk). (7) 

k 

The state |z/fc, p^) describes a neutrino with mass and momentum p^. It is a product of the 
states |z/fc) and \pk) belonging to the flavor- mass and momentum Hilbert spaces, respectively: 

Wk,Pk) = Wk) \Pk) ■ (8) 

In general the momenta of different mass eigenstates can be different, and they are deter- 
mined by energy-momentum conservation in the production process [^. For example, in 
the case of a muon neutrino produced by a pion decay at rest, 

7r+ ^ /i+ + z/^ , (9) 

the momentum of the k^^ mass eigenstate is 

4 V ml 2 V ml Ami 



^ Since detected neutrinos propagate along a macroscopic distance, we consider only one spatial dimension 
along the neutrino path. 
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with the corresponding energy given by 



^'-~V~^J +l"l^"^J + I^' ^^^^ 

where and are the masses of the pion and muon, respectively. The mass eigenstate 
li'kyPk) satisfies the energy eigenvalue equation 

'H\uk,Pk) = Ek\uk,Pk) , (12) 

where Ti. is the free neutrino Hamiltonian. 

In practice neutrinos with energy smaller than some fraction of MeV are not detectable, 
either because their energy is below threshold in charged-current processes, or because the 
cross section in elastic scattering processes, which increases with neutrino energy, is too 
small. Hence, if all the neutrino masses^ are much smaller than 1 MeV, all the detectable 
neutrinos are extremely relativistic and it is convenient to use a relativistic approximation 
for the momentum and energy of the mass eigenstates. 

At zeroth order in the relativistic approximation all neutrino masses are considered neg- 
ligible and all the mass eigenstates have the same energy E, equal to the modulus of their 
momentum: pk = Ek = E. The value of E is determined by energy-momentum conserva- 
tion in the production process. For example, from Eq. (0) one can see that in the pion 

decay process in Eq. we have E = ^ ^1 — ^ j ~ 30 MeV. Since the energy Ek and 

momentum pk of the k^^ mass eigenstate are related by the relativistic dispersion relation 

El=pl + ml, (13) 

the first order corrections to the equahties pk = Ek = E depend on the square of the neutrino 
mass and in general we can write 

2 

HI 

Pk^E-^^, (14) 

where ^ is a quantity determined by energy-momentum conservation in the production pro- 
cess. For example, from Eq. (|10D one can see that in pion decay ^ = | (^1 + — 0-8- Using 
the relativistic approximation of the energy-momentum dispersion relation in Eq. (0), 

Ek^Pk + ^, (15) 

from Eq. (|14|) we obtain that the energy Ek of the mass eigenstate neutrino z/^ at first order 
in the relativistic approximation is given by 

2 

7T? 

£;,^E + (i-0^- (16) 
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Direct searches of effective neutrino masses in lepton decays gave the upper hmits m^^ < 3 eV, < 



190 keV, rui,^ < 18.2 MeV (see Ref. |31 ). However, the sum of the masses of hght neutrinos is constrained to 



be smaller than about 24 eV by their contribution to the total energy density of the Universe (see Ref. [pl| 
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The wave function corresponding to the initial flavor state in Eq. (J^) is 

= M) = Y.u:k {x\pk)H) = 5^t/:,e^^'=(— w) , (17) 

k k 

where we have taken into account the initial coordinate xp of the production process. The 
time evolution of the wave function |^/'a(x)) is given by the Schrodinger equation 

^ ^ \€{^. t)) = n \^Plix, t)) = J2 Kk e^^'=^^-^^^ E, , (18) 

where we have used the energy eigenvalue equation ([1^). The solution of the evolution 
equation is 

\€M) = 5^f/:,e^P^(™)-^'=(*-*-) , (19) 

k 

where tp is the production time. At a time different from the production time tp and a 
coordinate different from the production coordinate xp, the state ^)) is a superposition 

of flavor states. Indeed, expressing the mass eigenstate {uk) in Eq. (|19D in terms of flavor 
eigenstates, 

l^^fe) = 5^t/„fc|z/„) (20) 

k 

(this relation is obtained by inverting Eq. (^) using the unitarity relation in Eq. (|^)), one 
obtains 

\€M) = J2 ( ^f/:,e'^'=(--^-)-^^^(*-*-) U,)j Wfs) . (21) 

13 \ k / 

This is a superposition of different flavor states, which reduces to one flavor only at the 
production point x = xp, t = tp, where the unitarity relation in Eq. (0) guarantees that 

|^^(X = Xp,t = tp)) = 

The term in parenthesis in Eq. ( pT]) is the amplitude of z/q, — >■ i/^ transitions, which is 
obtained by projecting = xp,,t = tp,)) at the detection point space and time, xp, and 

^D, on the flavor state |z/^): 

A.^,,(L,r) = (i.^|V.^(x = XD,t = tD)) =5^f/:,e^^'=^-^'^^[/^,, (22) 

k 

where L = xp,—xp is the source-detector distance and T = tp,—tp is the neutrino propagation 
time. 

Here it is to be emphasized that the detection process is described only by the flavor 
state 1 1//?), with no need to express the detected neutrino state as a superposition of mass 
eigenstates with different momenta. This is due to the principle of causality, which implies 
that different mass eigenstates with different momenta created in the production process 
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determine their own kinematics in the detection process where they are detected. The 
corresponding different momenta of the particles participating to the detection process do 
not spoil the coherence of the process if they are within the momentum uncertainty due to 
the localization of the detection process (through the uncertainty principle). 
The probability of vp transitions is 



leading to 



P.^^.JL,T) = \A.^^.JL,T)\\ (23) 



(L, T) = J2 \Uak\'\Upk\' + 2Re K.Up.U^.U;^ ^i{v.-P,)L-iiE,-E,)T ^24) 

k k>j 



It is clear that the transition probability at the detector depends only on the neutrino mixing 
matrix and on the momentum and energy differences of the mass eigenstates, which are due 
to their mass differences. Invoking now the relativistic approximations in Eqs. ([I^) and ([T6|) , 
we obtain 



p,„^,,(L,T) = ^if/„.rit/, 



/3k\ 



|2 



+ 2Re J2u:kUpkUajU;j exp 

k>j 



Am?. ^ ^ Am? ■ 



(25) 



with the mass-squared differences Am^^- = ml — m|. 

Equation (|25|) gives the oscillation probability as a function of the distance L between 
a source and a detector and the neutrino travel time T. However, in real experiments, the 
neutrino travel time is not measured, but the source-detector distance is known. Therefore, 
the neutrino travel time T must be expressed in terms of the distance L. To accomplish this 
task it is necessary to go beyond the plane wave approximation employed so far and remind 
that the propagating mass eigenstate neutrinos must be described by wave packets. Indeed, 
a plane wave has an infinite space-time extent, whereas neutrino production and detection 
are localized processes, that are appropriately described by the wave packet formalism. We 
associate each mass eigenstate, which has well defined kinematical properties, with a wave 
packet. The spatial width of the neutrino wave-packets is determined by the maximum 
between the quantum uncertainties of the position and time of the microscopic production 
process (see Ref. Q). The wave packets propagate with the group velocity 

z;. - ^ ~ 1 - ^ (26) 



where we have used the relativistic approximations (14) and (16). The different mass eigen- 
states can interfere if they overlap upon their arrival at the detector or if the coherence time 
of the detection process is larger than the time separation between the mass eigenstates 
0, 1^, 1^. In any case, from Eq. ( PB[ ) one can see that the arrival time at the detector of 
the mass eigenstate wave packets is equal to the source-detector distance^ plus a correction 

We use natural units in which c = h = 1. 
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proportional to m\/E'^. Since the phase in the neutrino oscillation probability in Eq. (|25|) is 
already proportional to Am1j/E, the correction to the relation T = L gives a contribution 
of higher order, which can be neglected. One may object that since this correction is also 
proportional to the distance L, it is not negligible when {Am\^L / E){m\/ E"^) > 1. This is 
true. However, one must observe that in practice neutrino oscillations are observable only 
when the phase {AmIjL/ E) is of order one. This is due to the fact that all neutrino beams 
have a spread in energy, all detectors have a finite energy resolution and the source-detector 
distance has a macroscopic uncertainty due to the sizes of the source and detector. When 
the oscillation phase is much larger than one, the averages over the neutrino energy and the 
source-detector distance wash out the oscillating term. If all the phases are large and all the 
oscillating terms are washed out, there remains only a constant transition probability given 
by the first term in Eq. 



{P.^^.,) = }^\UaknU^k\'. (27) 

k 

Therefore, the observability of neutrino oscillations guarantees that {Am\jL/ E){m\/ E"^) ^ 
1 and the correction to the relation T = L is negligible. 

Another correction to the relation T = L is due to the finite time interval in which 
the wave packets overlap with the detection process and to the finite coherence time of 
the detection process. One should average the oscillation probability in Eq. (p5|) over the 
resulting total time interval A^d- However, since the time interval A^d is characteristic of 
a microscopic process, it is always much smaller than the macroscopic propagation time 
T ~ L. If this property were not satisfied, oscillations could not be observed because 
they would be washed out by the average over time. As mentioned above, since neutrino 
oscillations are observable only when the phase {/S.m\-T / E) is of order one, the correction 
{/S.m\j/S.t-£)/ E) -C 1 is always neghgible. In other words, the phase of the oscillations is 
practically constant in the time interval At^ and the average over such a time interval is 
equivalent to the approximation T = L. 

From Eq. the appropriate approximation T = L leads to the transition probability 



P,^^,^{L) = J2\Uak\'\Upk\' + 2ReJ2u:kUpkUa,U;^ exp 

k k>j 



2E 



{2i 



as a function of the source-detector distance L measured in real experiments. This is the 
standard oscillation formula derived more than 20 years ago ^ Notice that as 



a result of the T = L relation, ^ has magically disappeared, leading to the independence 
of the neutrino oscillation probability from the specific details of the production process. 
From Eq. (^) one can see that there is an oscillating term associated with each pair of mass 
eigenstates k > j and the corresponding oscillation length is 

Since the state in Eq. (^) describing a flavor antineutrino has complex-conjugated el- 
ements of the mixing matrix with respect to the flavor neutrino state, the probability of 
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i^a '^13 oscillations is given by Eq. ( p8D with complex-conjugated elements of the mixing 
matrix. The oscillation lengths of neutrinos and antineutrinos are equal. 

Some comments are in order for the clarification of the derivation of the neutrino oscil- 
lation probability in Eq. (|28|): 

1) The initial neutrino state in Eq. is a coherent superposition of mass eigenstates. This 
is possible because neutrino masses are small and the differences between the momenta 
of different mass eigenstates is smaller than the momentum uncertainty of the production 
process due to its locahzation in space. Denoting the spatial interval of the source by 
Axp, we have the momentum uncertainty App ~ 1/Axp and the two mass eigenstates 
z/fc, Uj are emitted coherently if \pk — Pj\ ^ App. Using Eq. ( [T^ ) one gets the condition 
|Am^j | < 2EApp/^. For example, assuming Axp < 1 cm in the process in Eq. (P) of pion 
decay at rest, we obtain the coherence condition |Am^j| < 10^ eV^, which is well satisfied 
by neutrino masses below eV ranges. 

If the momentum uncertainty in the production process is too small, different mass eigen- 
states are emitted incoherently by the source and neutrinos do not oscillate. This happens 
if the momenta of the particles in the production process are measured with such a high 
accuracy to the extent that one can determine which neutrino mass eigenstate is pro- 
duced. In this case Axp > 2E /^\Amlj\ ~ Ll^j'^, i.e. the uncertainty in the localization 
of the source is larger than the oscillation length 0. It is obvious that in this case it 
is not possible to measure neutrino oscillations. Analogous considerations apply to the 
detection process. 

2) We have assumed that all the mass eigenstates are produced at the same space-time point 
(xp,tp) and detected at the same space-time point {xp),tp,). Since both the production 
and detection processes have uncertainties in space and time, different mass eigenstates 
can be produced and detected coherently at different locations in space and time within 
the uncertainties, generating interferences responsible for neutrino oscillations. However, 
one must pay special attention to the coherent character of the production and detection 
processes. In other words, there is a well-defined phase relation between the wave functions 
of neutrinos produced (detected) within the spatial and temporal coherence intervals of 
the production (detection) process. 

Let us first consider each mass eigenstate neutrino z/^ produced at a different space-time 
point (xp,tp), with |a;p — xp\ < Axp and |tp — tp\ < Atp. Since the component of the 
neutrino field that creates the corresponding mass eigenstate in the production process 
oscillates in space-time with a phase determined by the energy Ek and momentum pk, 
the initial phase difference between the wave function of the k}'^ mass eigenstate neutrino 
produced at (xp,tp) and that produced at (xp,tp) is 

A0P = Pk (4 - ^v) - Ek (4 - tp) . (30) 

Then, the phase at {x, t) of the wave function of the k^^ mass eigenstate neutrino produced 
at Xp and time tp is 

0fe(x, t) = A4 + Pk{x - x^) - Ekit - 4) =pkix- Xp) - Ekit - tp) , (31) 
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i.e. the same as that of the wave function of the k^^ mass eigenstate neutrino produced 
at xp and time tp, used in Eq. (pj|). 

Let us consider now each mass eigenstate neutrino z/^ detected at a corresponding space- 
time point {x'^,t^), with |xq — xd| ^ Aa;D and I^d — ^dI ^ AtD- The incoming neutrino 
wavef unction is 

A; 

Since the component of the neutrino field z/^ that destroys the corresponding mass eigen- 
state in the detection process oscillates in space-time with a phase determined by the 
energy and momentum p^, its phase difference between the points (x^, t^) and (xd, ^d) 
is 

A<^^ = Pk (x^ - 0:0) - Ek (t^ - to) . (33) 

This implies that the detected neutrino in a. oscillation experiment is no longer 

described by a simple flavor state |z/^) as in Eq. (0), but by the state 

= Y.Uke'''''-W) . (34) 

k 

which takes into account the coherence of the process. The amplitude of — > z^/3 transi- 
tions turns out to be 



which is equal to the standard amplitude in Eq. (P4 ). This mechanism allows the in- 
terference of different mass eigenstates even when the corresponding wave packets do no 
overlap, but arrive at the detector within the coherence time interval of the detection pro- 
cess 0, H, 1^. Clearly, in this case the different mass eigenstates are detected at different 
times, but the coherence of the detection process allows them to interfere. 

3) In the simple plane wave approach that we have adopted above, there is no indication 
of a coherence length for neutrino oscillations ||38[. The reason for the existence of a 
coherence length is that each massive neutrino propagates as a wave packet with its own 
group velocity. The overlap of different wave packets decreases with increasing distance 
from the source, until eventually they separate. For distances larger than the coherence 
length, two wave packets arrive at the detector separated by a time interval larger than the 
coherence time of the detection process. In this case the corresponding mass eigenstates 
cannot interfere and the oscillations are suppressed. Hence, in general there is a coherence 
length associated with each pair of mass eigenstates. At distances larger than the largest 
coherence length one can only measure the averaged constant probability in Eq. (p7|). This 
can be calculated using the wave packet formalism in the framework of quantum mechanics 
P, H or quantum field theory |l^]. However, one can estimate the coherence length 



in the following simple way, as done in Ref. [3S]. Using the group velocity formula in 
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Eq. (p6D, the separation of the k^^ and j^^ mass eigenstate wave packets after travehng a 
time T corresponding to a distance L is 

I Am?- 1 

\Axkj\ = \vk -Vj\Tc^ L . (36) 

The size of the mass eigenstate wave packets is given by the maximum between the 
spatial and temporal coherence widths of the production process. Denoting by the 
maximum between the size of the mass eigenstate wave packets and the spatial and 
temporal coherence widths of the detection process, the interference is suppressed for 
lAxfcjl > cr^., leading to the coherence length 

L™.^ ~ • (37) 

For n massive neutrinos there are n(n— 1)/2 oscillating terms in the probability in Eq. (|2 
each one with its oscillation length in Eq. (^) and coherence length in Eq. 



Several controversial and often misleading statements and derivations of the neutrino 
oscillation formula have recently appeared in the literature. Let us comment on some of 
them: 

A) The traditional derivation of the neutrino oscillation probability [jl| is based on the as- 
sumption that all the mass eigenstates have a common momentum pk = Pj = E. This 



is often called "equal momentum" treatment of neutrino oscillations. From Eq. (14) one 
can see that the equal momentum assumption is equivalent to assuming ^ = 0, i.e. this 
case is simply a particular case of the general kinematical relations in Eqs. (|1^) and ([T6|). 
Since the derivation of the neutrino oscillation formula in Eq. (pSD does not depend on 
the value of ^, it is clear that the equal momentum assumption is an acceptable one. 

B) It has been argued [El] that neutrinos created by a definite weak-interaction process must 



have a definite flavor at the production point at all times. From Eq. ([T9| ) one can see that 
this is possible only if all the mass eigenstates have the same energy, Ej. = Ej = E. Based 
on this considerations, the authors of Ref. claim that the "equal energy" assumption 



is the correct one. However, since the production process has uncertainties both in time 
and space, one could claim as well that neutrinos created by a definite weak-interaction 
process must have a definite flavor over all the spatial region of the production point at an 
initial time (or during the coherence time). From Eq. (jlUj) one can see that this is possible 
only if all the mass eigenstates have the same momentum, = pj = E, justifying the 
equal momentum assumption. Since energy and momentum satisfy the dispersion relation 
El = p\ + m\, different mass eigenstates cannot have simultaneously the same energy and 
momentum. Hence, it is clear that the equal energy assumption is as arbitrary as the equal 
momentum assumption, being yet another special case, ^ = 1, of the general kinematical 
relations in Eqs. ([I^) and (|l^). Obviously, the equal energy assumption leads to the same 



oscillation formula in Eq. (pH]) that we have derived without arbitrary assumptions. 
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C) It has been claimed that there is a "factor of two ambiguity" in the oscillation lengths 
p^ , |30| . This ambiguity arises from the observation that the propagation time T in 
Eq. (0) is not common among the different mass eigenstates, but one should take a 
different propagation time Tk for each mass eigenstate, given by Tk = L/vk, where Vk 
is the group velocity in Eq. (0). In this case the phase in the k-j interference term is 
4'kj = {Pk~Pj)L — {EkTk — EjTj). In the relativistic approximation E^Tk ~ pkL + miL/ E, 
leading to (pkj = '^^Ij^/^- C)ne can see that the phases are twice as large as the standard 
ones in Eq. (pS]), giving oscillation lengths that are a half of the standard ones in Eq. (p9|). 
This is the claimed factor of two ambiguity^. The confusion in this derivation is that it 
does not take into account the fact that different mass eigenstates arriving at different 
times at the detector must be absorbed coherently by the detection process in order to 
interfere. Hence, there is a phase relation of the detection process that must be taken 
into account, as discussed in item ^ above. This restores the missing factor of two in the 
oscillation lengths. 

D) As a variant of the equal momentum and equal energy assumptions, a "equal velocity" as- 
sumption has been proposed |2^, [Tl|, 0. Simple kinematics does not allow this possibility 



45[| . Indeed, Vk = vj implies that Pk/Ek = Pj/Ej. Taking the square of this relation and 
using the energy-momentum dispersion relation = El — m\, one finds mk/rrij = Ek/Ej. 
This equality cannot be true in a real experiment, since Ek/Ej ~ 1, whereas mk/rrij may 
be extremely small or large. 

E) It has been argued that a wave packet treatment of neutrino oscillations is unnecessary 
when neutrinos are emitted by a source that is stationary in time In this argument 
there is confusion between microscopic and macroscopic stationarity. In practice most 
neutrino sources (as the sun, or a reactor, or a supernova) can be considered to be sta- 
tionary sources of neutrinos. However, the microscopic processes that generate neutrinos 
in these sources cannot be considered stationary, because each microscopic process is able 
to emit neutrino waves only for a finite interval of time, the coherence time that we have 
denoted above by Atp. Hence, to understand the physics of neutrino oscillations, a wave 
packet treatment is necessary. We have seen above that a wave packet description of 
the propagating mass eigenstates is crucial for the determination of the correct relation 
between the neutrino propagation time and the source-detector distance. The necessity to 
establish such a connection is avoided in Ref. |^ by assuming that all mass eigenstates 
have the same energy. However, we have shown in item ^ above that in general such 
assumption is not justified (the stationarity argument given in Ref. obviously does 
not apply to the microscopic production process). 

F) In neutrino oscillation experiments neutrinos are usually produced through charged- 
current weak processes together with a charged lepton (see, for example, Eq. (|D). It 



'"^ A similar ambiguity has been claimed to exist in the case of kaon oscillations in Refs. [^9[ 
The authors of Refs. |4^, tried to solve this problem arguing that the interference between the Kl and 
Ks components should be calculated at the same time and the same space point. However, this argument 
does not hold because, as discussed in Ref. the coherent character of the detection process for a finite 
space-time interval allows wave functions at different space-times to interfere (see item above) . 
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is natural to ask if also the charged lepton has some kind of oscillatory behavior. In 
principle, flavor oscillations of charged leptons would be possible if their masses were very 
close so as to be produced coherently, and if there were some other means, besides mea- 
surement of mass, to distinguish different charged lepton flavors. However, since these two 
requirements are not satisfied in the real world, there are no flavor oscillations of charged 
leptons. Notice, in particular, that charged leptons are defined by their mass, which is the 
only property that distinguishes among charged lepton fiavorsQ. In other words, flavor is 
measured through mass. That is why charged leptons are defined as mass eigenstates. 



In Refs. [p6|, |28[| it has been claimed that the probability to detect the charged lepton 
oscillates in space-timej]. This claim has been refuted in Ref. 0]. The authors of Ref. p6| 
presented a proof that the oscillation in space-time of the probability to detect the charged 
lepton is a consequence of the oscillation in space-time of the probability to detect the 
neutrino. But it is well known that the probability to detect the neutrino, irrespective of 
its flavor, does not oscillate in space-time. This property is usually called "conservation 
of probability" or "unitarity" and is represented mathematically by the general relation 
Puc^fp = 1- Hence, the argument presented in Ref. actually proves that charged 
leptons do not oscillate! 



In conclusion, we have presented a simple but general treatment of neutrino oscillations 
in a quantum mechanical framework using plane waves and intuitive wave packets arguments 
when necessary. We have shown that the standard neutrino oscillation formula in Eq. ( pSD 
can be derived without arbitrary assumptions as "equal momentum" or "equal energy". 
We also tried to clarify some confusing statements on the quantum mechanics of neutrino 
oscillations that have recently appeared in the literature. 
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